clc;
clear;
[result,x]=Fi(100)
plot(x);